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Abstract: Gravitational waves were copiously produced in the early Universe whenever 
the processes taking place were sufficiently violent. The spectra of several of these gravita- 
tional wave backgrounds on subhorizon scales have been extensively studied in the litera- 
ture. In this paper we analyze the shape and amplitude of the gravitational wave spectrum 
on scales which are superhorizon at the time of production. Such gravitational waves are 
expected from the self ordering of randomly oriented scalar fields which can be present 
during a thermal phase transition or during preheating after hybrid inflation. We find 
that, if the gravitational wave source acts only during a small fraction of the Hubble time, 
the gravitational wave spectrum at frequencies lower than the expansion rate at the time 
of production behaves as Hgw(/) °^ f 3 with an amplitude much too small to be observable 
by gravitational wave observatories like LIGO, LISA or BBO. On the other hand, if the 
source is active for a much longer time, until a given mode which is initially superhorizon 
(fery* <C 1), enters the horizon, for kr) > 1, we find that the gravitational wave energy 
density is frequency independent, i.e. scale invariant. Moreover, its amplitude for a GUT 
scale scenario turns out to be within the range and sensitivity of BBO and marginally 
detectable by LIGO and LISA. This new gravitational wave background can compete with 
the one generated during inflation, and distinguishing both may require extra information. 
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1. Introduction 

Gravitational waves (GWs) are produced in the late Universe via cataclismic astrophysical 
events like hypernovae and inspiralling binaries. Because gravity is so weak, it is extremely 
difficult to detect directly with present day interferometers Q. On the other hand, dur- 
ing the violent processes which we expect took place in the very early Universe, several 
stochastic backgrounds of GWs of significant energy may be produced, although their am- 
plitude today is drastically reduced by the expansion of the Universe, making them equally 
difficult to detect (2|. Their discovery may however be possible in the near future, opening 
a completely new window into the uncharted territory of the very early Universe. For this 
we must determine the detailed GW spectrum, which strongly depends on the physical 
processes generating them. 

In the last few years there has been significant progress in the experimental prospects 
for detecting GWs with interferometers like LIGO and VIRGO and the future satellite 
mission LISA. This has stimulated research for sources of primordial GWs from the early 
Universe, either from hypothetical first order phase transitions [||, ||, [|, ||, [7], |8| or from the 
process of reheating after inflation g, [R], [TT], |T|, |T|, [jj, ||, |l§. 



The mechanism responsible for GW production during these early Universe phenomena 
is typically a causal process, like bubble collisions or turbulence, giving rise to spectra 
which peak at wavelengths that are well within the causal horizon during their generation. 
Thus, most of past analyses concentrate on contributions of GWs with wavelengths smaller 
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than the horizon at the time of production, with the exception of those generated during 
inflation 17], which are stretched by the inflationary expansion. 



In this paper we study the infrared behaviour of the GW spectrum produced either 
during preheating or during first order phase transitions, on scales which are superhorizon 
at the time of formation, i.e. k < Tt*, where k and are the comoving momentum 
and inverse horizon. We want to study a causal process of symmetry breaking like hybrid 
preheating 18, 19, 20, 22, 23], where the order parameter has global O(N) symmetry 



in the false vacuum and, upon symmetry breaking, the N fields undergo self-ordering on a 
given scale as soon as they enter the horizon, in particular on scales much larger than the 
inverse mass of the field in the true vacuum. 

We consider a multi-component scalar field which obtains a non-zero vacuum expecta- 
tion value (vev) v and a mass m, during a symmetry breaking process. We shall assume 
that this mass m is much larger than the Hubble parameter at the time of the tran- 
sition, since if the vev in the true vacuum is much smaller than the Planck scale, then 
H* ~ mv/Mp <C m. Such a model could describe the symmetry breaking process which 
triggers the end of hybrid inflation or a thermal phase transition. As long as we are only 
interested in superhorizon scales, k 3> TC*, we can neglect the radial, massive mode and 
treat the dynamics within the non-linear sigma-model (NLSM) approximation. On large 
scales, the anisotropic stresses are determined by gradient energy and the typical (comov- 
ing) scale is simply the time dependent horizon scale TL^ 1 . The field self-orders at the 
horizon scale, and the source of GWs decays inside the horizon. For scalar metric pertur- 
bations this process has been studied e.g. in Ref. [24]. It is very closely related to the 



scaling of global topological defects [25] even though for a number of components N > A 



there are no topological defects associated with such a scalar field in 3 + 1 dimensions. 

We work in the large N approximation within which the scalar field equation of motion, 
for scales larger than the inverse mass, k <C m, can be solved analytically. The GW 
spectrum will then be estimated by analytical approximations, introducing the anisotropic 
stress tensor sourced by the field fluctuations at different scales. 

Tensor perturbations from a NLSM in the large N approximation have also been stud- 
ied in Ref. p6| , p7] ], see also [28]. There the authors have calculated the tensor perturbation 



spectrum for scales which enter the horizon in the matter era and they have compared this 
with the inflationary signal in the CMB. Here we shall concentrate on the radiation dom- 
inated era and the detection of the signal in direct gravitational wave experiments like 



advanced LIGO pi, LISA pj and BBO pi]. 



The paper is organized as follows. In the next section we describe the formalism, derive 
the scalar field solutions and calculate the unequal time anisotropic stress correlators which 
source GWs. In Section || we study the production of GWs from long wavelength modes 
of this source. We derive a general formula that can be applied to different situations, 
depending how long the GW source is acting. In Section |4] we use this result to determine 
the shape and amplitude of the GW spectrum in two situations, first the case of a source 
producing GWs only during a small fraction of the Hubble time and, second, the case in 
which the source producing GWs acts for a much longer time, until a given mode which 
is initially superhorizon, krj* <C 1, enters the Hubble radius, krj ~ 1. In Section ||] we 
summarize our results and conclude. 

Notation Throughout this paper we assume a spatially flat Friedmann Universe with 
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metric 

ds 2 = a 2 (r]) (-d V 2 + SijdofdaP) , (1.1) 

where 77 denotes conformal time and we normalize the scale factor to unity today, 0(770) = 1. 
The comoving Hubble rate is TC = a' /a, while H = a' /a 2 is the physical one. The prime 
denotes derivative w.r.t. conformal time 77. 

2. Formalism 

We first introduce the NLSM and the large N limit of a global 0{N) symmetric scalar 
field, then we study the physics of the correlators of the anisotropic stress tensor. 

2.1 The model 

We consider an iV-component scalar field with a Lagrangian 

/ 2\ 2 

c = c Q + c l = -d^d^-\U T $- v -\ +A, (2.1) 

where $ T = 4>2, <I>n)/V2, A is the dimensionless self-coupling of <1> and v is the vev 
in the true vacuum. In the case of a thermal bath at high temperature, the Lagrangian 
Co obtains corrections of the form C\ ~ — T 2 <3? 2 , so that its minimum is at $ = which 
respects the global 0{N) symmetry of the Lagrangian. At low temperature, T < T c ~ v, 
the thermal corrections are too small to the keep the minimum at $ = and the global 
O(N) symmetry is spontaneously broken to 0(N — 1). In the context of hybrid preheating, 
there is no need for thermal restoration of the symmetry. The field $ acquires a large mass 
during inflation through its coupling to the inflaton x, £1 = ~9 2 ^ ^X 2 ■ Above a critical 
value, x > Xc = V^v/g, the effective quadratic mass of $ is positive and the field is fixed at 
$ = 0. When the quadratic mass becomes negative, x < Xc a tachyonic instability triggers 
the end of inflation and symmetry breaking. Soon after the symmetry is broken, thermal 
corrections and tachyonic effects can be neglected, and $ is closely confined (in most of 
space) to the vacuum manifold, given by ^ a cj) 2 (x, 77) = v 2 . Nevertheless, in positions 
such that their comoving distance is |x — x'| > H , the values $(x, 77) and Q(x.',rj) are 
uncorrelated, which leads to a gradient energy density associated to the iV — 1 Goldstone 
modes, p ~ (di&) 2 . For N > 2, the dynamics of the Goldstone modes is well described by 
a NLSM J32] , p5[ where we force ^ a 4> 2 = v 2 by a Lagrange multiplier. This corresponds 
to the limit A — > 00 in the above Lagrangian. This approximation is very good for physical 
scales with are much larger than m _1 = 1/(v / Atj). Of course, on small scales the field 
fluctuations still oscillates around the true vev, but in this paper we only focus on the 
superhorizon modes which are free to wander around in the vacuum manifold, giving rise 
to a gradient energy density which will generate GWs on these scales. 

Normalizing the symmetry breaking field to its vev, (3 = $/«, each component of the 



field obeys the non-linear sigma model evolution equation [24| 



Uj3 a - {dpP ■ d^l3)[3 a = , (2.2) 

where (8^(3 ■ d^fi) = ^ a r]^ d il f3 a {yi, rf)d lJ f3 a {^, 77) and ^ a /3 a (x, 7;)/3 a (x, 77) = 1. In the large 
AT- limit, we assume that the sum over components can be replaced by an ensemble average, 

T(x) = t V d^ a d v f3 a = N( V ^d^ a d u (3 a ) = f(n) . (2.3) 
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By dimensional considerations, T oc TL 2 , or 



f (77) = T oV ~ 2 , (2.4) 

with T D > 0. Replacing the non-linearity in the sigma-model by this expectation value we 
obtain a linear equation which can be solved exactly. In Fourier space it reads 

Pf + ~Pk + (k 2 ~ 5f) P a k = , (2.5) 

where 7 = d log a/d log 77 and primes denote derivatives w.r.t. 77. In a radiation dominated 
Universe 7 = 1 while in a matter dominated Universe 7 = 2. The solution to Eq. (2.5) for 
constant 7 is given by 

/? a (k,77) = (krfiv^foMMW + CkMYvikTi)] , (2.6) 

where 

1/2 =Q-7) +^0, (2.7) 

and Ci, C2 are constants of integration. Thus, v > 1/2 for a radiation dominated Universe 
and v > 3/2 for matter domination. Since in general we have that u > 0, Y u diverges for 
small argument, so we will keep only the regular mode of the solution J u , which can be 
written as 

1 



„ M) .Vl(^- T ^( M ,), (2. 



where f3 a (k, 77*) is the a-th component of the field at the initial time 77*. We assume that /3 is 
initially Gaussian distributed with a scale-invariant spectrum on large scales and vanishing 
power on small scales 



f (2vr) 3 C^5(k - k') , k V * « 1 



(/3 a (k,r7,)/3* b (k',^)) = I V> ^— Q u ^ " £ > 1 ( 2 -9) 

This means that the field is aligned on scales smaller than the comoving horizon 77* and 
has arbitrary orientation on scales larger than 77*. The condition that /3 2 = 1 actually 
introduces correlations between the different components of (3 but these lead to corrections 
of order 1/N to the above expression which we will neglect here. We also do not enter into 
the details of the decay of this function around kr]* = 1. The constant C is chosen such 
that the normalization condition is satisfied (up to corrections of order 1/N), 

I3 2 (x,r,,) = (/3 2 (x,^,)>(l+0(l/JV)) 

In the large TV- limit we neglect the corrections of order 1/iV which come from the fluctua- 
tions in f3 2 . On large scales this is a very good approximation. However, on small scales, 



- 4 - 



and in particular, on scales comparable with the inverse of the mass of the symmetry break- 
ing field, m _1 , the fluctuations are certainly not negligible. In our analysis we consider 
only large scales, where the above approximation is valid. 

In order for (/? 2 ) to be time independent we need that the equal time correlator be 
fixed to one: 

^* A {j) I dyy^4(y) = l, (2.11) 

where we have substituted C = Qir 2 r]l and we have set y = kr\. Note that the upper limit 
is actually ry/r/*, but at late times, the (dimensionless) integral is insensitive to the upper 
boundary, so we can take it to infinity and thus make the integral free of any time scale. 
In order to obtain a time-independent vev , we then just require 

i/ = 7 + l. (2.12) 
Introducing this relation into Eq. ( |2.7D , one obtains T a in terms of 7 as 

T = 3(7 + 1/4) . (2.13) 
The constant A is determined then by the condition 

l^jf^-^M, fcence A= V(»-W*-m . (2 , 4) 
Since v = 7 + 1, we can also write the amplitude of the field fluctuations, as 

n^ri) = Vl(£) 3/2 ^/3 a (M*) • (2-15) 
2.2 Unequal time correlators 

From Eqs. (^^) and ( |2.15| ) we obtain the following expression for the unequal time corre- 
lator of the field: 

7*bnJ J\\ - a ( VV'\ 3/2 Jv{kri)J v (k'i) / sn*bnj 



= (2 7 r)3 67 r^(W) 3/2 J ff^f; ^ - kO 

= (2vr) 3 5(k-k')^(A;,r / ,r / '). (2.16) 

We assume that the field (3 is Gaussian distributed initially. As its time evolution 
is linear, it will remain a Gaussian field and we can determine higher order correlators 
via Wick's theorem. This will be important in the next section when we determine the 
unequal time correlator of the anisotropic stresses which source the production of GWs. 
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Furthermore, this source is totally coherent [£E| in the sense that its unequal time correlator 
V% (k, 77, rj') is a product of a function of 77 and 77', 

V?M) = ^^^(W) 3/2 ^ff^y - , (2-17) 



with f(k, 77) = V6tt 2 A k 3/2 



Note the k 3 ! 2 scaling law at horizon crossing (krj ~ 1) which is characteristic for quantum 
fluctuations from de Sitter, i.e. inflation. This already hints to the fact that we will find 
a scale-invariant spectrum also in this case. 

3. The production of gravitational waves 

In this section we derive a general formula for the GW power spectrum sourced by super- 
horizon modes of a self ordering field. We also comment about the frequency range for the 
GW background produced in this way. 

Let us consider tensor perturbations (GWs) of the metric, 

ds 2 = a 2 (r/)(^ + 2h llv )dx»dx v , (3.1) 

where hij is traceless, h\ = 0, and divergence free, d' l hij = 0. Linearizing Einstein's 
equations yields the evolution equation of GWs sourced by the anisotropic stresses of the 
scalar fields <3>, 

hij (x, 77) + 2H h' t] (x, 77) - V 2 h tJ (x, 77) = 8vrG U {j (x, 77) , (3.2) 
where represents the TT part of the (effective) anisotropic stress tensor 

Ty(x, 77) = Wx, V)djr(*, V) ~ \&H [V<T(x, t/)] 2 . (3.3) 



Fourier transforming the GW evolution equation (3.2) we obtain 

h^(k, 7/) + 2H h' i3 (k, 77) + fcX(k, 77) = 8ttG A ijM (k)T lm (k, 77) (3.4) 
where the projector 

Pij(k) = Sj,j — kjkj , k = k/k , 
filters out the TT part of the Fourier transformed effective anisotropic stress tensor 

* f d 3 q 

n^(k,77) = A ijM (k) J ^3 qi q m <j) a (<i,r])(P a (k- q,7?) . (3.5) 

Note that we are summing over repeated indices both in coordinates and in field compo- 
nents. 
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The 2-point correlation function of the tensorial part of the anisotropic stress-tensor 
is of the form 



n^(k,r ? )n; m (k , ,r / ')) = (27r)^(k-k0n 2 (M^')^Wk) 



where 



Since the trace M. 



MijimQz) = - Ai jtlm (k) + A ijM (k) 



1. 



(3.6) 



(3.7) 



(3.f 



To determine H (k,rj,rj r ), we compute ^IIi J -(k,^)II^(k / ,77 / )y explicitly using Wick's theo- 
rem to reduce 4-point functions of the field to products of 2-point functions 



n^k^n^k',//)) = (2 7 r) 3 5(k-k')n 2 (A;,r / ,r ? ') . 



n ji (k,^)n? m (k' J f/) 

A-ij t pq(lt) A; mirs (k ) 



d?q d 3 q' 
(2tt) 3 (2tt) 3 



q P q q q' r q's{<t> a (<l,v)<i> a Q< ~ w)<Ftf,rf)<rQt ~ ^v' 



d 3 ad 3 



qaq , T 



q T Aq) ii (q' T Aq') im [^°(q,T 7 )0* a (q - ^(/(-qW^k' ~ + 



(2tt) 6 vj "" 

+ (^(q, ?? )^(q , ,r ? '))(<A a (k-q,r ? )0* 6 (k'-q',r ? ')) + 
+ (0 a (q, V W\k' - q', rf)) (0 a (k - q r/)^ b (q', r/))" 

(<? T Ag).. (q' T ^q') lm [v?(\<lU V)-P?W, rf, rf) <5(k)«5(k' 
+ (|q| , 7?, V)P?(|k - qU, 7?') *(q " q'M(k - q - k' + q') 
+ (|q| , 77, r?')7>f (|k - q|, 77, 77') 8{cf + q - k')5(q' + q - k) 



(3.9) 



where we use the notation (g T Ag) . . = qiAij i m q m and we have introduced the reality 
condition </>*(k) = 0(— k) and the unequal time correlator of the field </> which is defined in 
the same way as the one for /3, 



(^(k,77)0* 6 (k',77')) = (2vr) 3 «5(k-k')<(A;, 77, 77') . 



(3.10) 



The zero-mode of the anisotropic stresses vanishes due to isotropy so that the first term in 
the square bracket of the integral ( |3.9| ) does not contribute. 

We now can compute the unequal time correlator (Ily (k, 77)n* (k ; , 77') V Using 



( 9 T Ag)..( 9 T A g ).. = ^(l-(k.q)^ 2 



(3.11) 



we obtain 

U 2 (k, 77,77') 



d 3 



1 „4 



(2nf 



l-(k- q ) 2 l Vt(\q\,Ti,Ti')Vt(\k-q\,r,,r,'). (3.12) 
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We now relate the GW energy density spectrum to the unequal time anisotropic stress 
spectrum of the source, H 2 (k, rj, rf). For this we first write the GW evolution equation in 
momentum space, 

tiL + 2—h'u + k 2 hij = torGIUj . (3.13) 
J a J 

Defining a new variable hij = ahij, one obtains 

h'lj + (k 2 - ?-J hij = 8^Ga% . (3.14) 

In a radiation dominated background (a oc rj) this reduces to 

h'lj + k 2 hij = 8irGaUij ■ (3.15) 

The solution of this differential equation with the initial conditions hij = h\ - = is given 
by the convolution of the source with the Green function G(k,rj,rj') = sin(krj — krf), 

hij far) < rfcn) = —grj d V a(y/k)Uij(k,y/k)sm(x - y) , (3.16) 

where we have set x = kr/ and y = krj' . The source of gravity waves is acting for a time 
interval 5r]* = (r]n n — 77*) = e??*. If e < 1 we call the process short-lasting. This is the 
relevant case for example for GWs produced during a symmetry breaking phase transition 
where the source disappears after the phase transition since the latter typically lasts only 
for a fraction of the Hubble time. However, the Goldstone modes considered in this work 
may very well be long lived as they are not expected to interact with ordinary matter. In 
this case therefore a long lasting source may be better motivated. We discuss both cases 
below. 

After the source has decayed, GWs are freely propagating, and thus described by the 
homogeneous solution of Eq. fl3.15|) . 

hij(k, rj > r/ fin ) = A^(k) sin(fcr? - fcry fin ) + Bjj(k) cos(kr) - krj &n ) . (3.17) 

The coefficients Ay and Bij are fixed by matching the homogeneous solution to the inho- 
mogeneous one at rj = r]£ n . Matching both hij and its derivative h'^ yields 



8vrG 



1 in 



^i( k ) = —pr \ d y a (y/ k ) u ij(^ y/ k ) cos Ofin - v) 



8vrG f Xfin 

B ij( k ) = ~p~ J d V a {y/k) I hj(Ky/k)siD.{x &n - y) . (3.18) 



The GW energy density is given by (see e.g. ||) 

dpGVJ _ k 3 \h'\ 2 (k,r]) 
dlogk~ 2(2vr) 3 Ga 2 ' 

where the GW power spectrum has been normalized as follows: 



(3.19) 



h[j%ri)h'*j{ci,ri)) = 2(/ l ' + (k,r / )^(q,r ? ) + / l ' x (k,7 ? )/ l / ; :(q,r / )) = (2^) 3 5 3 (k-q) \h'\ 2 (k, V ) . 

(3.20) 
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Here our normalization differs from that of Ref. [27]. Their definition of the power spectrum 
is related to ours by 

V(k,7]) =2-Kk 3 \h\ 2 (k,rj) (3.21) 

and they infer = ~ whereas we obtain, with ( |3.19 ) and h' = kh for sub- 

horizon modes, 

dn GW (k, Vo ) _ k 5 \h\ 2 (k, v ) _ k 2 v(k, v ) 

dlogk 6ir 2 H 2 12ir 3 H 2 ' 

This difference in the normalization, which we attribute to an error in Ref. [^7]], leads to a 
reduction of the final result by about a factor 60, which may be relevant for observations. 
With the solution for hij above, we obtain for ij > r/fi n 

k 2 + n 2 /8ttg\ 2 r fln , r fln 



2a 2 V k 2 



dy J n dza (|) a (£) cos(z - y)Ii 2 (k, |, |) ,(3.22) 



where we have used Eq. (qJ^) . The GW energy density at time r\ is of course well defined 
only for waves with a wavelength well within the horizon, k S> H. Therefore we shall 
approximate k 2 + TL 2 ~ k 2 in the following. 

The GWs are sourced by the anisotropic stress of the scalar field 4> a = vf3 a . The 
correlators are simply related by 

Vf = v 2 Vf . 

With Eq. (|1|) we obtain the following expression for the GW energy density after the 
decay of the source, r\ > ?7fi n , 

dpcw(M) Gv 4 k 3 /"*» [*»» 

— ' or I dt, a{T)a(£) cos(k£ — kr) 



dlogk I- 1 a 1 !,/) ., 

x J d 3 pp 4 S m' l 6Vf(p,T,0Vf(\k-p\,T,0, (3.23) 

where cos 9 = kp. Inserting the power spectrum of (3 in the above expression and summing 
over the field components, we find 

d PGW (k,r]) Gv 4 k 3 36ir 4 A 2 [n* 



dlogk ' 4vr 4 a 4 (??) N 

x 



/ dr d£ a(r)a(£) cos(/c£ — kr) 



dWsm 4 ^ 3 ^ M ! ^ = , P l T) # = ■ (3-24) 
ik f<V^ (pr)" (PC)" (|k-p|r)" (Ik-pIO" 

Here the constant A comes from the normalization of (3, and it is given byEq. (pl^) . In 
the radiation dominated background considered here, we have v = 1+7 = 2 and A = 57r/4. 
Note also that we choose the normalization of the scale factor such that (1(770) = 1. Hence 
the comoving wave number k is simply related to the present frequency of the GW by 

/=-*-. 

J 2tt 
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In the next section we evaluate the present amplitude and frequency dependence of 
the GW spectrum generated in this way explicitly. For this, the following relation between 



temperature and time in a radiation dominated Universe are useful [34] 



Assuming an adiabatic expansion, g e g(aT) 3 = const., one finds 

, = -Mn- (aM) 1/3 L « V" _ i. 6 x io' S ec (™) >2*m ■ (^) 

T(r,)T„\ 2 J \,tf M (,)j V T ) y <* * ' v ' 

On the other hand, the expression for the temperature associated to a global 0{N) sym- 



metry breaking is [ 33 1 



/ 24 

T. = ^— 2 v , (3.27) 

independent of the coupling A. 

Before moving to the evaluation of Eq. ( [3.24 ), let us briefly determine the frequencies 



for the GW sources discussed in this paper. We are studying the IR modes k-q* < 1 of the 
GW spectrum, corresponding to frequencies smaller than the expansion rate at the time 
of production, /* = %*/(2tt), 

10 " (iv) Hz • < 3 - 28 > 

For the EW scale this corresponds to /,f w ~ 10~ 6 Hz, while for the GUT scale the asso- 
ciated frequency is f^ UT ~ 10 8 Hz. For a given energy scale M ~ T* at the time of pro- 
duction, we are describing one frequency range or another, but always frequencies smaller 
than the one corresponding today to that energy scale, / < f*(M) ~ 10~ 8 Hz(M/GeV). 
Clearly, only processes taking place in the radiation dominated Universe generate GWs 
with sufficiently high frequencies such that they can be observed by direct GW detection 
experiments. Indeed the frequency associated to the horizon at the matter-radiation equal- 
ity is far too small, /| q ~ 1CU 17 Hz, to be observed by direct GW detectors, like LIGO, 
LISA or BBO will be working. Therefore we consider only processes in the radiation 
dominated Universe and 7 = 1 and v = 2 are assumed for the rest of the paper. 

4. The gravitational wave spectrum today 

In this section we study two different cases, first the situation in which the source producing 
GWs lasts only a small fraction of the Hubble time at the moment of production and, 
second, the case in which the GW source acts for a much longer time, until the moment at 
which a given mode enters the horizon. 

4.1 Short lived source 

We first estimate the amplitude of the GW spectrum for large wavelengths, k < H*, from 
a short lived source which lasts from rj* to r/gn, such that (r/fl n — ?7*)/r/* = e <C 1 (as e.g. for 



the radial mode of <f) in hybrid preheating [12, |14|). Let us first note the following facts: 
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1) From Eq. ( [Oil) we see immediately that for small wavenumbers, kr]^ n <C 1, the re- 
sult scales like 

dlog 

2) Since the source is short lived, 77* « 77fi n , and we deal with superhorizon modes, fc?7* <C 1, 
we may set cos(kr] — krf) « 1 and the time integral can be replaced simply by a factor er;*. 

3) To estimate the momentum integral, we use that Bessel functions at small arguments, 
x = krj < 1, can be approximated by J v {x) ~ (x/2) u /T(l + i/). To obtain the dominant 
contribution at large wavelength (i.e. the least blue part) we may also set |k — qjT/* ~ qr)*. 

Using all the above considerations, we are left with a simple integral for the evaluation 
of the spectra of the IR modes (krj* <C 1) of GWs, at any time ij 3> 77* for which those 
modes have already crossed the horizon 



dpcwiv) 



dlog k 



Gv\ 4 2 k 3 2tt f 1 4 [V* / 

-367T A — / a cos # sin / dp 



t„«, It" « 4 (>7)JVJ-1 Jo '22T1(„+1) 

"^-W^ 2 = |^^(^) 4 ^(^) 3 , (4.1) 



II* 



where we used A = 5ir/4, v = 2 and we approximated J nfin dro(r)r 3 ~ a(r] ie )r^,8r] if = 
ea (v*) T li: since we have set r/gn — r?* = Srj* ~ er/*. 

With this we can now evaluate the ratio of the GW energy density to the critical 
density today, for the IR modes kr]* <C 1, as 

1 dpGw(r) ) 5vr 4 / v \ 4 e 2 3 

°GW(/) = TT— 7— « TT^8 77- ^rad(^) 



p c dlog A; 7- 2 8 V^W N 



~ 10 ' s {jtSi (kn - )3 ' (4 - 2) 

where we used Hi = 8ttGp^/3, we expressed the radiation density today as p ra( j ~ 
p*(a*/ao) 4 and we introduced the the radiation density parameter today as $l ra d ~ 4.2 x 
1CP 5 . We have also neglected the factors coming from the ratio of the effective relativistic 
degrees of freedom since they appear only with the power 1/3. 

Note that this formula is general for the IR spectrum of GWs generated at any process 
in which the source, a iV-component scalar field, has rapidly acquired its true vev v at 77* 
and undergoes a short phase of self-ordering which lasts for a fraction e < 1 of the Hubble 
time. 

Finally, note also that very generically we have 77* oc T^ 1 oc 1/v so that f^cw <x 
A: 3 oc v fc 3 and not as v 4 , as one could naively have concluded from Eq. 



4.1.1 The electroweak phase transition 

The comoving horizon size at the electroweak (EW) phase transition is given by the EW 
energy scale T* ~ 100 GeV, g c g(T*) = 106.75, 

7/* ~ 7.5 x 10 sec . 
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Inserting this above with / = fc/ (2vr), we find 



GWU; 7 -2 s Iad \M F iJ A^V mHz y N\mRz) V 7 

For the last expression we have used v ~ T*. This result is of course unmeasurably small. 

4.1.2 A GUT scale phase transition 

To have any chance to measure this spectrum, we need a vev which is not too many orders 
of magnitude below that Planck scale, since the GW energy density is suppressed by a 
fourth power of the ratio of the vev to Mp\. The best change might be a GUT scale with a 
vev of the order of v ~ 10 16 GeV. But then of course 77* will be very small and the dominant 
contribution will come from very high frequencies, lower frequencies being suppressed by 
the factor (fer/*) 3 . For T* = 10 16 GeV we have 

??* ~ 5 x 1CT 10 sec , 

leading to 

Apart from the fact that this result suffers severe additional suppression at measurable 
frequencies which are significantly below 1GHz = 10 9 Hz, the sensitivity of 10 _12 il rac j ~ 
10 -16 cannot be reached with any presently proposed experiment at those frequencies. 



Therefore, we can only conclude that the superhorizon GW spectrum generated from 
a short lived self ordering scalar field is much below presently proposed experimental sen- 
sitivities. 

4.2 A long lived source 

As we have seen in the previous subsection, short lived Goldstone modes cannot lead to a 
significant GW background. But since Goldstone modes are typically non- interacting and 
long lived, it is more natural to consider them for a time which is much longer than the 
horizon scale r/*. To compute the GW energy density produced by such a self ordering 
scalar field, we consider Eq. ( [3.24 ) and set r]n n = rj k = since the solution ( |2.15| ) decays 



inside the horizon, when kr] > 1. We then have to compute the following integral 



dp GW {k,r] k ) _ Gv 4 k 3 36tt 4 A 2 
dlogk ~ l^a 4 ^.) N 



rl/k rl/k 

/ dr I d£ a(r)a(£) cos — kr) x 



/ 



p £^ | < 1 <1 ^ (prY (pty (|k-p|r)" (Ik-pIO" 



Note that the range of integration of the variable p in the above expression is set to be 
{pr}* < 1, |p — k|r/* < 1} since the initial two point correlator of the scalar field turns out 
to be different from zero only in this range of momenta [c.f. Eq. (2.9)]. 



In order to obtain an analytical result for the above integral, we perform the following 
approximations : 
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• We are interested in scales k that are superhorizon for all the time of GW production, 
namely kr < 1 and fc£ < 1 for times r, £ between 77* and r/fin = l//c, therefore we 
approximate cos(/c£ — kr) ~ 1 . 

• We neglect the angular dependence of |p — k| so that the angular integral reduces to 
2vr / sin 4 6»(icos(9 = 32vr/15. 

• In the range of integration where pr > 1 we substitute |k — p|r ~ pr, while when 
pr -C 1 we approximate |k — p|r <C 1. 

• The range of momenta for which we can expand the Bessel functions in terms of small 
arguments isp < min(l/r, l/£), while in the range min(l/r, l/£) < p < max(l/r, l/£) 
we should distinguish between large and small argument expansions of the Bessel 
functions. Finally, in the range max(l/r, l/£) < p < 1 / 77* one can consider the large 
argument limit for all the four Bessel functions of the above integral. 

Taking into account all the above considerations, we find that the complete integral 
becomes 



/ dr (% dpf(p,r,0 

Jn* Jfj* JO 




which allows us to separate the integral in p using the asymptotic behaviour of the Bessel 
functions, 

x u 

J v (x) ~ ; for I « 1 , 

, , Pi ( (2u + l)vr\ 

Jv\x) — \ — cos x for x 3> 1 . 

V xtt y 4 J 

We can distinguish three different intervals: 

• The IR contribution, I±(k), for < p < 1/r, with |k — p|r < 1 and |k — p|£ < 1 . 

• The mixed (UV+IR) contribution, l2(k), for 1/r < p < l/£, with |k — p|r ~ pr > 1 
but |k — p|£ ~ p£ < 1 . 

• The UV contribution, ^(/c), for l/£ < p < I/77*, with |k — p|t ~ pr > 1 and 
|k — p|£ ~ p£ > 1 . 



Therefore we can finally write 

dpGw(k,r] k ) 
dlog A: 



V{k)[h{k) + I 2 {k)+h(k)] , (4.6) 



where the pre-factor T>{k) contains the coefficients in front of the integral in Eq. ([4.5|) , the 
factor coming from the angular integration (32tt/15) and the factor 2 that comes from the 
symmetry of the double time integration, namely 

Gv 4 k 3 36vr 4 A 2 32vr n Gv 4 k 3 1C , 3 

V(k) = — -, — j—— x x 2 = — -15 • 4vr 3 . (4.7) 

1 ; 4vr 4 a 4 ( Vk ) N 15 iV a 4 r/ fc 1 ; 
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The three integrals of Eq. (|4.6| ) are given by 



h(k) 



l/k 



4096 
4096 A; 3 35 





it)* 




nl/k 


J, 


(It 






1 


ri 



-3 t 3 



1A 



dp p 



-4 tf 4 



1/r 



dp p 



3 - g(^*) 3 + 2( fcr ?*) 5 



d,(pr) J„(pg) J„(|k - p|r) J„(|k - pJQ 
(*"-)" (pf)" (|k-p|r)" (Ik-pIO" 



(4.8) 



h{k) 



and 
hik) 



rl/k pt 

/ dr d£ o(r) a(£) r 3 £ 3 



dp p f 



32vr 



dr I di r 4 e 4 



l/T 

l/r {PT) 



,6 



6 J„(pr) J„(pQ J„(|k - p|r) J„(|k - p|0 
(pr)" (p£f (|k-p|«r)" (|k-p|£)" 
5tt" 



dp p - > 
5 cos ( pr 



128vr k 3 



(4.9) 



l/k 



dr I di o(t) a(£) t 3 £ 3 



4ffp^rad 
7T 2 

3tt 2 /fc 3 



l/fc 



dr / d£ r 4 £ 4 



dp p' 



6 J,(pr) J,(p£) J„(|k - p|r) J v (|k - p|0 



(pt)" (pO u (|k-p|r)" (|k-p|£) 1 



dpp 6 2 / 5vr\ 2 / 5vr 
2 cos I pr — ] cos ( p4 — 



/c (p?-) 5 (p6 



4 / 



1 



1 



9 9 



(/cr/*) 3 - (£;?7*) E 



-log 2 ^*) - -log(A;77*; 



4 

(4.10) 



More precisely, in the above computation we substituted each cos 2 x by its mean value 
(cos 2 x) = 1/2 averaged over a few oscillations, and we introduced the usual expression 
for the scale factor in a radiation dominated background, a(rj) ~ Hq^Qj-^t], which is 
consistent with a$ = I today. 

All three terms have a scale-invariant spectrum. Actually, the "UV" contribution 
given in Eq. ( [4.10 ) is the largest. Summing all the three contribution and considering 
the dominant part in the limit fcry* <C 1 [hence also (fer/*) 3 log(/cry*) <C 1], we obtain the 
following scale-invariant spectrum 



dpGw(k,f]k) 
dlog k 



~ 5 • 2 V 



^rad Pc 



1 



+ 



1 



+ 



1 



2 12 • 105 2 6 tt • 45 27vr 2 



60 x 



^rad Pc 



(4.11) 



where we have used the Friedmann equation Hq = 8irGp c /3. Redshifting the above ex- 
pression until today, we obtain for the GW energy density parameter, 



^Gw(Mo) 



dpGw(k,rj ) dp GW (k,rj k ) 4 „ 

"« \Vk) — — iZ rad 



Pad log k 



p c dlogk 



60 

N 



(4.12) 
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Figure 1: The sensitivity of present and future GW experiments are compared with our results for 
a long lasting source and inflation. We show, the amplitude of the scale-invariant GW background 
expected from a GUT scale inflation (blue, dashed) and from a self-ordering long lived source 
as studied in this paper, for a symmetry breaking field with N = 4 real components and a vev 
v = 10 -2 Mpi (top, red line), v = 10 _3 Mpi (middle, blue line, overlying with inflation) and v = 
10~ 4 A/pi (bottom, green line). The big dot at the right end of the horizontal lines represents the 
frequency ( 3.28| ) associated to the horizon at the initial time of production. 



This corresponds to a scale-invariant GW spectrum produced by a self-ordering scalar field 
in the large TV- limit. This result is valid for all wave numbers k which enter the horizon 
when the Goldstone modes of our iV-component field are still massless and the field has not 
yet decayed. Scales which enter the horizon after this time rjfa, i.e. scales with krj^ n < 1, 
are suppressed by a factor (kr]^ n ) 3 , as for them the result for a short lived source with 77* 
replaced by r/fi n applies. 

4.3 Numerical integration 

In order to obtain more accurate results, and to check the validity of our analytical ap- 



proximations, we have also performed a numerical evaluation of the integrals in Eq. ( 3.24 ). 
If we set the final time of integration to be the horizon crossing, r]^ n = l/k, as we did in 
the analytical evaluation for the long lasting source (f4.5|), we obtain the following result 
for the final GW density parameter today 

^Gw(Mo) - ^^rad (~^~) ' ??fin = 1 A • (4.13) 

This suggests that the analytical approximation somewhat overestimates the result. How- 
ever, we can continue the integration to later times when the wavelength has already 
entered the horizon. 
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Figure 2: The density parameter in gravitational waves as a function of krj. For scales outside the 
horizon, kr] < ir, we observe the (kr]) 3 dependence (short dashed line), while for scales that have 
entered the horizon, kr) > it, the GW energy density saturates, at a normalized value of 511 (long 
dashed line). This result implies a significant scale-invariant GW spectrum today. 



The integral in Eq. ( |4.5| ) allows us to compute the GW energy density in the limit 
<C 1, using the 

GVft rad 



Na A (r]) 



7517*1 dqq z F(q) 



kr) 

dx cosx J 2 (qx) 



kr, 




1 2 


" rkr) 


+ 


\ dx 




Jo 



dx sin x J 2 (qx) 



where the kernel F(q) comes from the integration over angles, 



2\2 



1 



du (1 - u 1 ) 
(q 2 + 1 - 2qu) 2 ~ 2V 



16g + I2q{q 2 - I) 2 + 3(q 2 - lf{q 2 + 1) log 



(4.14) 
(q + l) 2 



and we have made the approximation, J2{x\/q 2 + 1 — 2qu) — > J2(qx), inside the time 
integration. We have checked that for large times the result is correct within 0.1%. 

Numerically evaluating ( [4.14 ), we find that the GW energy density continues to grow 
until horizon crossing, krj ~ ir, and saturates thereafter, see Fig. H This agrees with the 
result of Ref. who find a peak in the power spectrum V(k,rj) at approximately this 
value, and also explains the l/a(rj) 2 dependence of the Power spectrum, V oc VLq^ /a 2 , for 
scales that have already entered the horizon. 

For krj S> 4 the gravitational wave energy density saturates at a value 

511 / v \ 4 

^Gw(Mo) - ^^rad ( , (4.15) 

where we used again the usual normalization of the scale factor in a radiation dominated 
background. These results suggest that the GW spectrum produced by this mechanism 
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still grows inside the horizon and reaches its final value somewhat after horizon crossing. 
This is consistent with the fact that the power of the scalar field that sources these GWs is 
not absent inside the horizon, but it is indeed given by the Bessel functions in Eq. ( 2.17| ), 
which decay rather slowly as functions of krj. 

In the following analysis we will consider the numbers arising from the numerical 
integration, as given in Eq. (4.15). 



4.4 Observational constraints 



Our result for the amplitude of the GW spectrum (4.15) is inside the range of detectability 



of the BBO H experiment (0 G w(&) > 10~ 17 ) and is marginally detectable by LISA |30|] 
or advanced LIGO @ (^Gw(k) > 1(T 10 ). Indeed, with ft rad ~ 4.2 x 10~ 5 , we find that 



BBO would detect this signal if the symmetry breaking scale v satisfies 

4 

>4.7-10~ 16 iV => > 1.5 • 10~ 4 iV 1/4 



Mpi J ~ m pi 
Concerning the sensitivity of LIGO or LISA, the signal is detectable if 

4 

> 4.7 • 10~ 9 N => > 0.008 N 1/4 . 



M Pl J ~ Mpi 

In other words, for scales higher or around the GUT scale, v > 10 16 GeV, the very long 
wavelength tail which we have studied here could be observed. 

In order to relate the above scale-invariant GW energy density to the GW spectrum 



from inflation, we compute the relative tensor-to-scalar ratio r. Following Ref. [35], one 
has the following expression for the GW density parameter from inflation 

" O w(fc,>»)=4.36xl0-r(£)" T , r s gg> , (4.16) 

where k = 0.002 /iMpc" 1 , V T (k) = rVs(k )(k/k ) nT and we used the WMAP result, 
'Ps(ko) = 2.21 x 10~ 9 . This concerns only the wavelengths which enter the horizon in the 
radiation dominated era, before equality. Comparing the above expression for ut — with 
our Eq. ( 4.15| ), we obtain in our case 



Another usefull comparison with inflation is the relative strength of the GW energy den- 
sities produced by the above two different mechanisms. Considering always wavelengths 
which enter the horizon in the radiation dominated epoch, we have p£ 



«- 1O - 13 (lO^) 2 - 8 - 4 x 1O - 5 (0- <« 8 » 

where M denotes the energy scale of inflation, H% = 87rGM 4 /3. The ratio between the 
GW energy density produced by our mechanism and the one from inflation is then 

«(Mo) n[m) 
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Comparing these results with those of Ref. [27], where the authors mainly concentrate 
on the spectrum of GWs produced in a matter dominated universe, we reproduce perfectly 
the amplitude of their spectrum V(k,r)) defined as in Eq. flOU) , but their final relat lve 
strength 1Z is nearly 2 orders of magnitude larger than what we find in Eq. ( [Hp . We 
believe this is due to the factor l/(2-7r 3 ) missing in their expression for HgwI^Io) which 
has to be introduced for consistency with the definition of the power spectrum V{k). 



5. Conclusions 



In this paper we have estimated the contributions to the gravitational wave background 
from a symmetry breaking phase transition on large scales, kr]* < 1. We have concentrated 
on the analysis of the Goldstone modes and we obtained the following main conclusions. 

If the modes are short lived with duration er/*, e < 1 their contribution is blue and 
suppressed by a factor e 2 (/c?7*) 3 . This result is actually generic, independent of the nature 
of the short lived source. Indeed, one typically obtains 

n G w(fc) - (k v ,) 3 n rad n 2 x e 2 , (5.1) 

where £lx is the density parameter of the source of anisotropic stresses at the moment of 
creation. For the Goldstone modes the factor £l 2 x is replaced by (v/Mpi) 4 . This strong 
suppression factor renders GWs from short-lived Goldstone modes entirely unobservable. 

The situation is different for long lived Goldstone modes. There the suppression factor 
(fery*) 3 is absent. Therefore, if the Goldstone modes remain massless until a time ?7fi n , for 
modes with krj^ n > 1 the spectrum is scale invariant and the amplitude is given by 

511 / v \ 4 

n G w(k)^ w n iaA (—J , (5.2) 

which is marginally detectable with the experimental sensitivity of advanced LIGO or LISA 
and is well within the range of BBO for a GUT scale phase transition. The results for the 
long-lived source are summarized in Fig. |l[ 

If the Goldstone modes are still present at decoupling, r/fi n > -q^ec , these GWs will also 
leave a signature in the cosmic microwave background where they lead to a scale-invariant 



contribution very similar to the one of global textures, i.e. a N = 4 global O(N) model [25]. 

Note that this new GW background from self-ordering fields after inflation (e.g. from 
hybrid preheating) has a power spectrum very similar to that coming from inflation, and 
therefore it may become important to disentangle both if they are present simultaneously, 
that is if the scale of inflation and that of symmetry breaking are related by parameters of 
order one, like in hybrid inflation. 
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